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ABSTRACT 

Frequency  parameter  plane  techniques  are  developed  for 
the  analysis  and  synthesis  of  both  linear  and  nonlinear 
systems.   The  transmission  function  of  a  dynamic  system 
contains  two  parameters  which  are  adjusted  to  provide 
desired  performance.   These  parameters  represent  elements 
of  the  system. 

The  four  dimensional  paremeter  space  relating  magnitude, 
frequency,  and  the  two  parameters  may  be  represented  by  seven 
two  dimensional  projections  of  the  parametric  curves. 
Computer  programs  are  developed  for  the  computation,  and 
graphs  of  results  are  obtained.   Interpretations  of  the  graphs 
can  be  used  for  design  of  systems  such  as  tank  stabilizers. 
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I.   INTRODUCTION 

The  analysis  and  synthesis  of  control  systems  may  be 
accomplished  using  the  frequency  parameter  plane  technique. 
A  primary  consideration  in  control  system  design  is  the 
adjustment  of  the  system  frequency  response.   To  achieve 
this,  the  frequency  parameter  plane  technique  may  be  employed 
and  is  based  on  a  graphical  procedure  for  considering  four 
parameters,  alpha,  beta,  frequency  and  magnitude  of  the 
transmission  functions.   The  objective  is  to  choose  values 
for  alpha  and  beta  (or  functional  relationships  if  the  design 
is  to  be  nonlinear)  to  provide  the  desired  frequency  response 
characteristics.   This  choice  is  possible  through  proper  use 
of  the  frequency  parameter  plane  curve  families. 

The  goal  of  this  thesis  is  design  using  analysis  techniques 
The  parameter  frequency  response  method  is  used  to  synthesis 
and  design  systems.   A  steam-engine  governor  is  used  in  the 
analysis  as  an  illustration.   The  antiroll  tank  stabilizer  of 
a  ship  is  also  designed  using  the  frequency  parameter  plane, 
for  both  passive  and  active  tanks.   The  tank  size  and  active 
device  characteristics  are  represented  by  the  parameters, 
which  are  varied.   Parameter  plane  curve  families  are  used 
in  the  design  procedure  presented  in  later  chapters. 


10 


II.   PARAMETER  PLANE  ANALYSIS 

Dynamic  systems  may  be  either  linear  or  nonlinear. 
Parameter  plane  analysis  permits  us  to  study  dynamic  behavior 
by  treating  system  parameters  as  variables.   One  first 
determines  the  characteristic  polynomial  of  the  system  which 
may  be  in  the  form 

k  N  N-l 

F(s)    =   E      ak   S      =   aNS'VaN_1SIM    '  + +   a^+ag    =   0 

K  ™  U 

(2.1) 
where   ak  =   coefficients   which    are    functions    of   the   parameters. 

S      =  c+jw' =    -cw   +jw      A      2~  =   u   ej9 

By  substituting  the  value  of  S  into  equation  2.1, 
equation  2.2  is  obtained. 

F(s)  =  E   a.  J<  ejke  (2.2) 

k=0   K   n 

Since   equation    2.2    has    complex   coefficients    it   can    be   converted 
into    two   equations    by    requiring    reals    and   imaginaries    to   be 
zero   independently. 

Note   that   Sk    =   ojk   ejke   =   wk    (Cos    ke   +   j    Sin    ke) 

=    Xk   +   J    Yk 
where    XR    =   wkCos    ke    =   wkCos(k    Cos"1?)    =    (-l)kwkTkU) 

Yk    =  u>kSin   ke    =   oJkSin(k    Cos-1?)    =    (-l)kojkUkU) 
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T|((^)  and  Uk(?)  are  Chebyshev  polynomials  with  values 
as  shown  in  Table  1,  2,  and  3. 


It  follows  that 
N 


k  k. 


Z  a.  [(-l)VTk(t)]  =  0 
k=0  K         K 


N 


k  k, 


Z   a-Tt-UVMc')]  =  0 
k=0  K        K 


(2.3) 


are  simultaneous  equations  which  may  be  solved  for  two 
undefined  parameters  in  the  a.  . 

For  many  practical  problems  the  algebraic  form  of  the 
coefficients  is  one  of  the  following  two  cases: 

1 .  ak  =  bka  +  ck3  +  dk 

2.  ak  =  bka  +  ck3  +  hkae  +  dk 

Case  1  can  be  solved  by  Cramer's  rule  and  the  results 
for  a  and  6  are 


a  = 


3   = 


C1D2  -    C2D1 

B1C2  _    B2C1 

B2D1  -    D2B1 

B1C2  "    B2C1 


where 


Bl    \lQ<-V Vn  V^       B2  -kfiot^Vn \W 


Cl    iM^Vn^k^ 


N  k         k 

D,    =    Z    (-D    dkco   kT.(?) 
i      k=Q  k    n      k 


N  k  k 

C?    =    S    (-D    C.03    KU.(?) 
c      k=[)  k    n      K 

°2  "J  «-1>kdk»nk,Jk.^ 
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Case  2  by  substituting  equation  2.3  rearranges  into 
the  form. 

aB1  +  6C1  +  agH1  +  D]  =  0 

aB2  +  BC2  +  a0H2  +  D2  =  0     ....(2.4) 

where  the  definitions  of  case  1  still  apply,  and 

H,  =jo(-l>VnkTk(?);  H2  =j0(-1'khk"nkuk'5) 

Eliminating   3   the   solutions   are   obtained 

_   -e±/ez-   4ac  t 0   c    % 

al,2 21 (2-5a) 

and   corresponding    values    for   3    are 

Blal    2  +  Dl  B2al    2+Dl 

<h>2  ■  -h,.,;^,  =  -  H2a,;2+c2    <2-5b> 

Eliminating   a,    the   following   solution   for  B 

_   -f±/f -4bd  (2   6   x 

Bl,2  2b  U'5a' 

and   the    corresponding    values    for   a   are 
_,     C,Bl  .g*Pl    _        C2Bl-^D2 
"1.2  Hl6l>2+Bl  H26li2+B2      <2-6b> 

where   the    notation    is    the    following 

a    =   B2H1    -    B.jH2;  c    =    C-,D2    -    C^ 

b    =    C2H]    -    C^;  d    =    B^    -    B^ 

e    =   B2C]    -    B1C2+H1D2-H2D1 
f    =    C2B]    -    B2C1+H1D2-H2D1 
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As  can  be  seen,  there  are  generally  two  pairs  of  solu- 
tions (a, 3)  in  terms  of  u  and  £.   The  solutions  are  given 
by  equations  2.5a  and  2.5b  provided  a  f   0 .   If  a  =  0  the 
solutions  are  given  by  equations  2.6a  and  2.6b  provided 
b  f   0.   When  a  =  b  =  0  case  2  reduces  to  case  1. 

Illustrati  ons  I .   A  feedback  control  system  is  as  shown 


R 


"v  o 

Kl 

C 

J     L 

S(S+l)(S+2) 

K-S 

2 

' 

The    closed-loop    transfer   function    is 

K. 


T(s)    = 


1 


S3   +    3S2    +    (2+K1K2)S    +    K] 


The  characteristic  equation  of  the  system  is 
S3  +  3S2  +  (2+K]K2)S  +  K,  =  0 


Let 


a  =  K^2 


6  =  K 


1 


The  equation  becomes: 


S3  +  3S2  +  (2+a)S  +3=0 


From  the  parameter  plane  method 


Bl  -■*>„ 

Cl  ■  1 


B2  =  -% 
c2  =  0 
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D1    =    (4C3-3?)(-wn3)    +    (2C2-l)(3a)n2)    +    (-2Cu>n) 


D2    =    (4?2-l)(-a)n3)    +    (20(3o)n2)    +    (-2u)n) 


SO,    the   solutions    for   a  and   3    are 

a  =— 4oj      C     +    6ai    C  +   (jo        -    2 

n    s  n^         n 

6   =    3^2    -   2?»n3 


The   c-curves    and   and   w   -curves    on    the   a-3    plane    are   shown 
on    Fi  gure   2.1. 

II 1 ustrati  on    II .      The   system   is    the   same    as    illustration    I 
let        K]    =    3 
K2        a 

The    characteristic   equation    becomes: 

S3   +    3S2    +    (2+a3)S    +3=0 
Then 


B]    =   0 

Cl    =    I 
Hl    =    "^n 


B2  =  0 
c2  =  0 
H2    =    -con 


D]    =    (3S-4.C3.)    ajn3   +    (6?"-3)wn2    -    2?wn 
D2    =    (1-4S2)    un3   +    6Cwn2    -    2a)n 


And   the    other   definitions    are 


a   =   0; 


b    =    «on; 


c   =   0 


d    =    -4oj    3C3      +  6    A2   +    (3(j3-2a3n)  X.   -    3u2 
n  n  n        n'  n 

e   =    f   =   2cj4  K  -  3w3 
n  n 
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So   only   one    solution    for  a   and   $    in   equations    2.6a   and   2.6b 
is   obtained: 

0   -    3.2    -   2a£ 


a  = 


*tf 


6co   C 
n^ 


2    .    0 
«h   +   2 


2w3C 
n^ 


3w. 


9r  Un     -     2 

2g   +   _n 

0 


The    curves    are   shown   on    Figure   2.2 

Consider  the  curves  on  Figures  2.1  and  2.2:  the  roots  of 
the  characteristic  equation  are  defined  as  functions  of  a  and 
3.      If  we    choose   3    =   40    and   a  =    20    on    Figure    2.1,    then 

C      =   0.111 
therefore,    the    roots    of   the   system   are: 

S    =    -   0.5    ±    J4.47 

and   the    corresponding    value    for   the    real    root   is    r2.      In    like 

manner   on    Figure    2.2,    if  we    desire    roots    at   x,   ~   0.111    and 

u      =   4.5   we    found   a   =    0.5    and   3    =   40. 

On   the   other  hand,    if  c   =   0.2   and  co     =   5   the    roots    are 

n 

-l±j4.9  and  -1.   From  the  Figure  we  find: 
3  =  K1  =  25 

a  =  25  (Figure  2.1)  or  1.0  (Figure  2.2)  therefore  K2  =  1  .,0 
Parameter  plane  methods  are  yery   convenient.   The  a-3 
plane  can  be  used  to  find  the  system  parameters,  and/or  roots 
of  the  characteristic  equation  at  any  defined  values  of  the 
parameter.   The  stability  of  the  system  can  be  investigated 
at  the  line  £=0  which  corresponds  to  the  imaginary  axis  on 
the  S-plane. 
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Table   2.      Functions   TkU) 
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Table  3.   Functions  U. (?) 
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Figure  2.1   Root  Parameter  Plane  Curves  for  illustration  1 
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Figure  2.2  Root  Parameter  Plane  Curves  for  Illustration  2 
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III.   PARAMETER  FREQUENCY  RESPONSE  TECHNIQUES 

A.  INTRODUCTION 

Both  linear  and  non  linear  dynamic  systems  may  be  sub- 
jected to  periodic  inputs.   It  is  important  therefore  that 
we  be  able  to  analyze  such  systems  in  the  frequency  domain. 
If  the  system  is  linear  there  also  exists  a  correlation 
between  frequency  response  and  time  response  as  guaranteed 
by  the  Fourier  integral  which  permits  us  to  estimate  the 
transient  response  from  the  frequency  response.   Thus  the 
frequency  domain  may  be  used  to  design  systems  for  non- 
periodic  inputs.   In  this  chapter  frequency  response  analysis 
of  linear  systems  on  a  frequency  parameter  plane  is  discussed 
i n  detai 1 . 

B.  ANALYSIS 

A  linear  system  can  be  described  mathematically  by  a 
transmission  function: 


T(s) 


_  N(s) 


Din 

where    the    system  may   be    either   open-loop   or   closed-loop   and 

N(s)    and   D(s)    are   polynomials.      In   polynomial    notation 

m  k 

I      a.SK 

T(s)    -   k=0   ' 


?      b    Sq 
q  =  0    q 


where,  a's  and  b's  are  coefficients  which  are  functions  of 
the  system  parameters. 
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If  we  choose  two  systems  parameters  a  and  $  and  if  we 
consider  only  sinusoidal  signals  the  frequency  response 
transmission  function  becomes 

T(JW)  -  D(jo).a.B) 

0  m 

aQ  +  a<|j'U  +  a2(jo>)  + +  am(jw) 

bQ  +  b^u  +  b2(ju)2+ +  bn(ja))n 

[aQ  -  a2w  +  a.u  -  ]+j[a-|U)-a3to  +  a5a)  -...] 

[bQ  -  b?w  +  b.u   - ]+j [b,o)-b_o)  +brW  -...] 

In  abbreviated  notation 

Ne(u,a,B)    +    jN0(co,a,B) 
T(ja))    =    Dglu.a.B)    +   jDo(a),a,B) 

Simi  larly 

Ne((i),a,B)    -    j    N0(a3,a,B) 

T("ju)    =    De(Usa,B)    -    j    Do(W,a,B) 
Let   the   magnitude   be    represented   by   M: 
M   =    |T(jo))| 

M2    =    |T(jo))|2      =       |T(Joi)|     |T(-Ju)| 

therefore 

2         [Ne(w,a,B)]2    -    [N    (u>,a,3)]2 

M      =   s ~ 

[De(ai,a,3)]^    -    [D0(w,a,3)] 

from  whi  ch 

2          4                   2                          3  2 

«        [a^-a^to  +a.w    - ]      +    [a-,w-a-co  + ] 

~                      2           4                     2                             3  2 

[bQ-bpW   +b.u    - ]      +    [b-iW-b-o)   + ] 

2m  T    •  '   ir        2m_1  "i-    '      t 

2  (-D    a.    w     +   Z  (-1)    *    a.     w 

_    k  =  0  (even)  K k  =  1    (odd)  K 

2n  g_  2n-l  gjj_ 

Z  (-1)2K'       q    +    E  (-1)       2       '      q 

q  =  0  (even)      bq   w  q=l    (odd)  bq   w 
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Then 


m 


M 


I      (-1)  at  u 
2   k=0 


k  "   2k 

k 


i   (-i)V  »2' 

q=0       H 


where,  the  a.  's  and  b  's  are  functions  of  a  and  g.   So,  there 
are  four  variables  M,  w,  a,  and  B  in  this  equation.   The 
equation  can  be  plotted  by  fixing  one  variable  and  using  one 
variable  as  a  parameter.   Seven  different  combinations  are 
possible . 

Exampl e .   Consider  the  system  of  a  Steam  engine  Governor 
[9]  for  which  the  transmission  function  can  be  written  as: 


C1C2 

I  m 


c       S4+(Cg+Ce)S3+(Ke+Kg+CeCg)S2+(KeCg+CeKg+ClC2)S+KeKg 
ml  Imlm  Imlmlm  Im 


where    <f>   =    instantaneous    speed    response. 
<f>     =   command   speed. 

m   =   equivalent   mass    of   the   governor   sleeve. 
C      =   damping    coefficient   at    the    governor   sleeve. 
K     =   stiffness    of   governor   spring. 
C,    =    increase    in    the    upward    force    on    the   governor 

sleeve . 
C2    =    increase    in    steam   torque. 
I    =   equivalent   moment   of   inertia   of   the   engine. 
C      =   damping   coefficient    of   the    engine. 
K      =   stiffness    of   engine    spring. 
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From  studying  the  mechanical  characteristics  [10], 
possible  values  of  the  parameters  of  the  system  are  sub- 
stituted in  the  equation.   Let  K  and  C2  be  varied  by  adjusting 
the  governor  spring  and  the  valves  for  steam  flow  respectively. 
Then  the  equation  becomes: 


4> 


6  26aS 


*c        S4+10.5S3+(30+25B)S2+(137.5B+6  26a+12.5)S+62.5B 
where        a  =    C„ 


3   =    K, 


then 


M2   = 


[626aa)]' 


[w4-(30+253)w2+62.56p+[-10.5a)3-(137.53+626a+12.5)a)]2 


From  the  curves  of  Figure  3.1  several  useful  bits  of 
information  can  be  obtained  by  inspection: 

(a)  At  any  intersection  between  two  constant  w  curves 
the  coordinates  define  the  a-3  pair,  and  the  values  of  u  of 
the  intersecting  curves  define  the  bandwidth  (since  the  mag- 
nitude is  -3  db)  . 

(b)  A  locus  of  constant  co  define  a  set  of  a,  B  pairs 
which  are  potentially  useful  in  self  adaptive  applications, 
i.e.,  if,  in  a  given  system,  B  is  subject  to  change,  values 
of  a  can  be  determined  such  that  a  magnitude  M  =  M,  can  be 
maintained  at  id  =  u,  ,  by  determining  the  values  assumed  by  6 
and  providing  a  suitable  adjustment  in  a. 

(c)  If  a  number  of  families  of  curves  are  generated  with 

coordinates  as  in  Figure  3.1  but  with  different  M  values,  the 

entire  frequency  response  can  be  determined  for  any  chosen 

a,  B  pai  r. 
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From  the  curves  of  Figure  3.2a,  b  and  c,  it  can  be  seen 
that: 

(a)  Each  frequency  has  only  one  a,  0  pair  which  gives 
the  maximum  magnitude. 

i.e.,  at  &j  =  10,  a  =  1.0,  8  =  3.0  (Figure  3.2a) 
u  =  12,  a  =  1.4,  3  =  4.5  (Figure  3.2b) 
a)  =  14,   a  =  1.75,  6  =  6.6  (Figure  3.2c) 

(b)  The  locus  of  the  maximum  magnitude  point  is  a 
straight  line  when  frequency  is  increased. 

(c)  The  axis  which  is  drawn  through  the  ellipses  gives 
the  values  of  6  (and/or  a)  and  magnitude  which  are  coordinates 
of  points  on  the  envelope  of  the  a  versus  magnitude  curve  as 
in  Figure  3.3. 

On  Figure  3.3  are  curves  of  a  versus  magnitude.   These 
curves  as  drawn  show  the  envelope  of  the  family  of  curves 
that  corresponds  to  the  axis  of  the  ellipses  on  Figure  3.2. 
This  envelope  bounds  all  possible  magnitudes  at  a  frequency 
of  12  rad/sec,  for  all  Q. 

Figure  3.4,  which  presents  a  family  of  constant  magnitude 
curves  on  a  versus  frequency  coordinates  (B  =  constant)  pro- 
vides a  global  view  of  the  possible  frequency  responses  of 
the  system: 

(a)  We  can  plot  the  Bode  diagram  by  reading  the  frequency 
and  magnitude  at  the  intersections  along  any  line  of  constant  a, 

(b)  The  maximum  magnitude  for  the  given  values  of  B  is 
found  at  the  center  of  the  family  of  ellipses,  and  also 
increases  as  a  straight  line  when  6  is  increased. 
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(c)   At  constant  3  (3=5.0),  the  curve  of  magnitude 
is  -3  db.,  it  shows  the  bandwidth  is  between  the  intersection 
of  the  curve  and  any  line  of  constant  a,  such  as  at  a  =  2 
the  bandwidth  is  10.5  rad/sec. 

A  plot  of  3  versus  frequency  with  constant  a  may  be  used 
in  the  same  way. 

Curves  of  a  versus  magnitude  are  shown  on  Figure  3.5 
with  3  constant  and  frequency  the  family  parameter. 

(a)  The  Bode  diagram  can  be  plotted  at  any  a  by  evaluating 
magnitude  and  frequency  along  a  horizontal  line. 

(b)  At  the  intersections  of  two  frequency  curves  the 
Bode  diagram  for  that  value  of  a  will  have  the  same  magnitude 
at  each  of  the  frequencies. 

(c)  Each  frequency  curve  defines  one  value  of  a  (3  is 
constant)  for  which  the  magnitude  is  maximum  at  that  frequency. 

(d)  If  two  frequency  curves  intersect  each  other  at 
the  magnitude  -3  db .  ,  the  baridwidth  at  that  value  of  a  is 
defined  by  the  two  frequencies. 

Figure  3.6  plots  a  versus  frequency  with  3  as  the  para- 
meter.  Each  3-curve  is  intersected  twice  by  each  a  line, 
thus  defining  two  values  of  frequency  for  the  one  value  of 
the  a,  3  pair,  these  two  values  of  frequency  give  the  band- 
width of  the  system  (since  the  magnitude  is  -3  db). 

Figure  3.7  shows  a  family  of  Bode  diagrams  with  a  as  the 
family  parameter  and  3  constant.  It  may  be  used  to  select  a 
value  of  a,  or  to  define  a  desired  variation  in  a  as  a  function 
of  frequency. 
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Each  curve  family  is  seen  to  contain  useful  information. 
Depnding  on  the  type  of  information  needed  some  curve  families 
may  be  more  useful  than  others. 

C.    STABILITY  ON  THE  FREQUENCY  PARAMETER  PLANE 

Figure  3.8  shows  the  stability  on  the  root  parameter  plane 
of  the  above  system.   The  curve  of  z,   =   0  is  obtained  from: 

1  .248  -  0.4902O32  +  5.541u3  +  0.19968oo4 
a  =  


u>(25w2  -  62.5) 


4    on    2 

a  _   a)   -  30u 


25a)2  -  62.5 


In  the  typical  system,  the  region  of  stability  on  the 
frequency  parameter  plane  can  be  determined  by  using  the 
above  relations. 

The  stability  boundary  (?=0)  curve  has  been  added  on 
Figure  3.1,  through  Figure  3.6  and  the  shaded  areas  are  the 
stable  region.   Synthesis  and  design  with  the  frequency  para- 
meter plane  are  considered  only  for  stable  systems. 

Each  Figure  has  a  different  interpretation.   Figure  3.1 
plots  a  versus  3  with  magnitude  constant  at  -3  db ,  and  frequency 
as  the  family  parameter.   By  choosing  an  operating  point  in  the 
stable  region  at  the  intersection  of  two  frequency  curves,  one 
guarantees  both  stability  and  bandwidth.   Figure  3.2  uses 
frequency  as  the  constant  parameter  and  magnitude  as  the 
family  parameter.   This  permits  optimum  choice  of  a  and  3 
for  desire  magnitude  and  frequency. 
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The  other  curves  can  be  considered  in  the  same  way.   At 

any  point  on  the  Figures  in  the  stable  region,  the  parameters 

at  that  point  can  use  to  design  the  system  and  obtain  values 
of  magnitude  and  frequency  as  desired. 
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Figure  3.1 


Frequency  Parameter  Plane  with  Stability  Limits 
Plot  of  a  vs  3  with  M  constant  at  -3  db  and  w 
as  family  parameter. 
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Figure  3.2a 


Frequency  Parameter  Plane  with  Stability  Limits. 
Plot  of  a  vs  3  with  w  constant  at  10  rad/sec 
and  M  as  family  parameter. 
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Figure  3.2b   Frequency  Parameter  Plane  with  Stability  Limits 
Plot  of  a  vs  6  with  to  constant  at  12  rad/sec 
and  M  as  family  parameter. 
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Figure  3.2c   Frequency  Parameter  Plane  with  Stability  Limits 
Plot  of  a  vs  3  with  a>  constant  at  14  rad/sec 
and  M  as  family  parameter. 
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Figure  3.3 


Frequency  Parameter  Plane  with  Stability  Limits 
Plot  of  a  vs  M  with  co  constant  at  12  rad/sec 
and  3  as  family  parameter. 
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Figure  3.4a   Frequency  Parameter  Plane  with  Stability  Limits 
Plot  of  a  vs  to  with  3  constant  at  3  =  4  and  M  as 
family  parameter. 
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Figure  3.4b   Frequency  Parameter  Plane  with  Stability  Limits 
Plot  of  a  vs  w  with  3  constant  at  3=5  and  M 
as  family  parameter 
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Figure  3.4c  Frequency  Parameter  Plane  with  Stability  Limits 
Plot  of  a  vs  co  with  £  constant  at  B  =  6  and  M 
as  family  parameter. 
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Figure  3.5  Frequency  Parameter  Plane  with  Stability  Limits 
Plot  of  a  vs  M  with  3  constant  at  6  =  5  at  co  as 
family  parameter. 
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Figure  3.6  Frequency  Parameter  Plane  with  Stability  Limits 
Plot  of  a  vs  03  with  M  constant  at  -3  db  and  3 
as  the  family  parameter. 
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Figure  3.7 


Frequency  Parameter  Plane  (Bode  Diagram).   Plot 
of  M  vs  u  with  3  constant  at  3=5  and  a  as  family 
parameter. 
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IV.      NONLINEAR   SYSTEM   TECHNIQUES 

Generally,   nonlinear  systems   consist  of  linear  parts 
which   are   described  by   linear  differential    equations    and 
nonlinear  parts  which   are   described  by   nonlinear  functions. 
If  the  system  is   nonlinear  with   an   applied  sinusoidal    input, 
the  signal    at   various   points   in  the   system  may  or  may   not  be 
sinusoidal.      In   general    the   signals  will    be  periodic,   but 
may   contain  sub-harmonics    and/or  harmonics. 

The  most  widely   used  method  in   the   stability  analysis 
of  nonlinear  control    systems   linearized  by   the    describing 
function   technique   is    the    Nyquist   diagram.      With   the    concept 
of  the   variable   critical    point,   which    depends    upon   the    frequency 
and  signal    level    at  the    input   of  nonli nearf ty ,    the   Nyquist 
diagram  permits   a   simple   stability   analysis    and  investigation 
of  sustained  oscillation. 

However,   the    Nyquist   diagram  is   inconvenient   to   apply   to 
multiloop  control    systems.      Difficulties   also   arise  with   the 
adjustment   of  the   system  parameter,   since   the    diagram  should 
be   plotted  each    time   a   parameter  is    changed..     Thus,    this 
method   can  hardly  be   applied   to   design   problems    in  which   the 
describing    function   does    not   appear  as    a   variable   gain   factor. 
Other  methods   may   be    used,    namely   the  ;Hurwitz    criterion, 
Root-locus    techniques,   Bode   diagrams,   Nichols    charts, 
Mitrovic's   method,    approximation  method,    Popov   method,    and 
phase   plane    techniques,    etc.    which    can   be    used   to   investigate 
the   nonlinearity   of  the   system.      However   the    parameter   plane 
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method  is  more  convenient  since  the  nonlinearity  values  can 
be  represented  by  the  parameters. 

Consider  the  sinusoidal  system; 

T(ja))"  D(ju>,a,S) 


Then 


N(ja),a,B) 
D(jw,a,3) 


M    =  |T(ju)|  = 
Define  one  parameter  as  B  to  be  constant  and  the  other  para- 


meter a  represents  the  value  of  the  nonlinear  element.   So; 

N(  jh),a) 


M 


D(jw,a) 

The  equation  has  three  variables:  magnitude,  frequency,  and 
the  nonlinear  parameter  a,  which  the  parameter  plane  method 

is  to  plot. 

The  block  diagram  of  the  nonlinear  system  may  be  the 

followi  ng : 

N 


■^ 

a 

G(s) 

J 

o — o 


G(s) 


a 


a  is  the  nonlinear  element. 

6(s)   represented  the  linear  system 
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Nonlinear  elements  such  as  square  law  devices,  coulomb 
friction,  dead  space,  saturation,  backlash,  relays,  etc.  may 
be  represented  by  their  describing  function.   The  parameter 
is  then  used  to  represent  any  adjustable  parameter  in  the 
describing  function,  and  the  resulting  plots  on  the  frequency 
parameter  plane  indicate  the  effect  of  the  nonlinearity  on 
the  system  behavior. 
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V.   DYNAMIC  SYSTEM  DESIGN 

A.    INTRODUCTION 

In  the  design  of  compensators  for  linear  or  nonlinear 
systems,  linear  system  design  to  make  the  system  stable  uses 
the  root  locus,  Bede  plot,  Nyquist  and  Nichol  chart.   However 
the  parameter  plane  method  is  more  convenient  for  both  linear 
and  nonlinear  systems.   In  classical  methods  of  the  root 
locus  and  frequency  domain  (Bode  diagram),  the  compensators 
are  introduced  to  reshape  the  corresponding  graphs  of  the 
original  system  in  order  to  achieve  desired  system  response 
characteristics . 

In  the  parameter  plane  approach  the  direct  effects  of 
the  additional  blocks  and  their  connection  "is  to  change  the 
characteristic  equation,  so  that  the  additional  adjustable 
parameters  are  conveniently  distributed  among  the  coefficients 
of  the  characteristic  equation.   This  enables  us  to  have 
control  over  the  characteristic  roots  by  varying  the  introduced 
parameters . 

The  parameter  plane  is  used  here  to  investigate  and 
design  the  system  for  antiroll  stabilization  of  a  ship  by 
means  of  passive  tanks.   The  size  and  properties  of  the  tank 
are  designed.   The  problem  is  extended  by  including  the  pump 
system  thus  converting  the  system  to  an  active  stabilization 
system. 
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B.    THE  DYNAMICS  OF  A  SHIP  IN  ROLL  MODE 

A  mathematical  description  of  the  dynamic  motion  of  a 
ship  requires  six  degrees  freedom,  three  equations  of  force 
and  three  moment  equations,  these  six  equations  are  too 
complex  for  analytical  solution.   If  only  the  rolling  of 
the  ship  is  considered,  it  may  be  modeled  as  a  pure  rota- 
tional motion  about  the  longitudinal  axis  of  ship,  with  no 
moments  about  the  Y  and  Z  axes. 

The  ship  in  rolling,  is  analogous  to  a  pendulum.   The 
dynamic  equation  of  a  pendulum  is: 

m  f  +  B  i  +  ^  <f>  =  F 
where  m  =  mass  of  the  pendulum. 
B  =  damping  coefficient. 
g  =  gravity  of  the  earth. 
1  =  the  length  of  the  pendulum. 
<J>  =  the  instantaneous  angle  of  pendulum. 
F  =  the  excite  force. 
Including  hydrodynamic  damping,  the  dynamic  rolling 
equation  of  a  ship  is  the  same  as  the  pendulum. 

s     y  s     y  s     y  .yi 

where      J      =   total    moment   of   inertia   around   axis    of   roll. 
B      =   hydrodynamic    damping    coefficient. 
K     =   static    righting   moment    coefficient. 
Y   =    angle    of   drift,    radians. 
K     =    drift-angle    coupling    coefficient. 
<f>   =    instantaneous    angle   of   rolling. 
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The  equation  can  be  written  as: 


r    ^s   s  T    sT    y 


B. 


where  c  =  damping  factor  of  the  ship.  =  /i — jt- 


w 


c    =    natural    frequency   of   the    ship    roll.    =V~r~ 
>  s 


When    the   ship   is    operating    in   a    sea   state,    there   is    a     > 
multitude   of   forcing    frequencies    present   and   the    integrated 
response    to   all    these    frequencies   must   be   considered.      Further 
experiments    are    therefore    carried   out   with    the   model    at   zero 
speed    in    irregular   beam  waves.      The    rolling   of   the   ship    is 
caused   by    the    force   due    to   the   wave   slope.      Then    the    equation 
becomes  : 

$   +    2    £s    ojs    <j>    +   ws    <J>    =   a)      if) 

where  $=    the  magnitude  of  wave  slope. 

C.    MODEL  OF  A  SHIP 

The  model  of  ship  used  for  the  example,  to  simulate  and 
design  the  tank  stabilizer  is  one  given  in  Reference  [6],  the 
data  and  characteristics  of  the  ship  were  the  following: 

Displacement  weight  w  =  936  tons. 

Metacentric  height  GM  =  0.73  meter. 


Ship  moment  of  inertia 


J   =  1290  tons-m-sec' 


Hydrodynamic  damping  coefficient    B   =  60  tons-m-sec. 
Static  righting  moment  coefficient  K  =  684  tons-m. 
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The  mass  moment  of  inertia  of  the  tank  stabilizer  (J^) 
increased  the  value  of  J  ,  but  it  is  very  small  and  can  be 
considered  negligible  when  compared  to  J  ,  for  convenience 
use  the  value  of  J   as  above. 

The  beam  of  the  ship  can  be  found  by, 


B    = 


2tt/GTT 

"V5 


2tt/5M 
Cto_ 


where      C   =   a   coefficient   dependent   upon   the   moment   of   roll 

and   approximate    is   0.797   sec/m2. 
Therefore    the    beam   of   the    ship   can   be   calcuiHated   and    the 
result   is: 

B   =   9.25  meters. 


D.    PASSIVE  TANK  STABILIZATION 
1 .    Analysis 

A  passive  tank  can  only  be  designed  for  optimum 
performance  for  one  ship  condition  of  displacement  and  meta- 
centric height.   When  the  ship  condition  varies,  the  free 
surface  correction  y  varies,  together  with  ifche  quantity  of 
water  required.   It  is  therefore  necessary  to  check  the 
performance  of  the  stabilizer  under  different  ship  loading 
conditions  and  it  could  be  that  a  compromise  design  has  to  be 
adopted  which  is  acceptable  for  all  conditions. 

The  free  surface  loss  is  generally  expressed  as: 

1  =   5  w  B 
V   48  A  h 
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where 


and 


i  =  the  second  moment  of  inertia  of  the  free  surface 

V  =  the  volume  of  displacement  of  the  ship. 

w  =  the  weight  of  active  fluid. 

A  =  displacement  of  the  ship. 

B  =  beam  of  the  ship. 

h  =  depth  of  fluid  in  tank. 

—  =  usual  notation  for  free  surface  correction, 
v" 


The  tank  size  parameter  is  limited  primarily  by 
static  stability  requirements.   However,  the  tank  should  be 
designed  with  as  large  a  free  surface  as  possible  consistent 
with  the  requirements  of  static  stability.   Choosing  values 
of  y  is  equivalent  to  varying  the  tank  size. 

Since  the  pendulum  system  is  analogous  to  ship 
rolling,  the  effect  on  the  ship  of  a  hydrodynamic  tank 
stabilizer  is  analogous  to  the  addition  of  a  second  pendulum, 
The  combination  of  two  pendulums  in  series,  may  be  called  a 
double-pendulum  system  as  shown  in  Figure  5.1. 


Figure  5.1   Double-Pendulum  System. 
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The  first  pendulum  represents  the  ship  rolling 
system,  and  the  second  pendulum  represents  the  dynamic  tank 
system.   The  length  of  the  pendulums  should  be  the  same,  so 
that  the  periods  are  nearly  equal  which  means  that  both 
natural  frequencies  of  the  ship  rolling  and  the  tank  are 
approximately  the  same.   When  the  ship  mass  has  been  swung, 
the  free  oscillations  are  observed.   When  the  tank  mass  has 
been  attached,  the  masses  of  both  ship  and  tank  can  be  swung 
and  the  interplay  and  exchange  of  energy  between  the  tank 
and  the  ship  can  be  seen. 

The  tank  stabilization  system  uses  water  or  some 
fluid  mass  device.   There  are  many  types  of  tanks.   They  can 
•be  designed  as  Doughnut  tanks,  completely  filled  sea  ducted 
tanks,  partially  filled  sea  ducted  tanks,  IE-tube  tanks  and 
free  surface  tanks.   These  tanks  systems  were  suggested  by 
Dr.  J.  H.  Chadwick  Jr.   Typical  tanks  such  as  free  surface 
and  U-tube  tanks  have  been  suggested  by  G.  J.  Goodrich  and 
Williams  C.  Webster  respectively.   Both  types  are  described 
and  designed  in  this  chapter. 

A  block  diagram  of  ship  and  tank  system  is: 
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For  a  passive  tank  system  the  "Control  block" 
does  not  exist,  or  may  be  considered  a  unit  gain.   Thus  the 
block  diagram  for  a  passive  system  becomes: 


* 


a 


SHIP 


* 


TANK 


The   relationship   of   the    ship   with    a   tank   stabilizer 
can   be    represented   as    a    double   pendulum.      The   equations    are: 


J*  +   Bc$     +   K<J>     +   J.e      +   Kc.e      =    F 


stv 


'st' 


Jst<(>  Kst<j>  +  Jte  +  Bte  +  Kte  =  o 


where     J    =  the  moment  of  inertia  of  the  ship  rolling. 
B    =  hydrodynamic  damping  coefficient  of  ship. 
K    =  the  static  righting  moment  coefficient  of 

ship. 
J.    =  the  moment  of  inertia  of  the  tank. 
B.    =  the  damping  coefficient  of  the  tank. 
K.    =  the  static  equivalent  coefficient  of  the 

tank . 
J  .   =  the  coupling  moment  of  inertia  of  ship/tank. 
K  .   =  the  coupling  static  coefficient  of  ship/tank. 
F-,    =  the  force  insert  in  ship  as  the  wave  slope. 
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2.    Free  Surface  Tank  Design 

The  free  surface  tank  as  Figure  5.2  shown,  is 
advantageous  in  decreasing  the  magnitude  gain  of  rolling, 
and  also  decreases  stability  of  the  ship  since  the  free 
surface  loss  is  high. 


Figure  5.2  Free  Surface  Tank. 

The  equations  of  the  ship  and  tank  can  be  rewritten 
as:  for  the  ship 

••  •      p        j   p         0  0 

<j>   +    2cswsct>   +   ol>s 4>  +   y-  tos9    -    ycos    =    w^ 
for   the    tank 


H        0  **  0  **  *  0 

£■  u)£<j)  -  w£<j>  +  e  +  2ctute  +  w^e  =  o 


where  w  =  natural  frequency  of  the  ship  roll 


to 


.  =  natural  frequency  of  the  tank 


C   =  damping  factor  of  the  ship  roll 
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?.  =  damping  factor  of  the  tank. 

2 
g  =  gravity  of  the  earth  =  9.8  m/sec. 

d  =  the  distance  between  the  water  level  in  the 

tank  and  the  roll  center. 

ty     =  wave  slope  function. 

y  =  the  free  surface  correction  factor. 
_L    _!_!>£ 

V  GM    2AGM 

Both  equations  have  the  negative  sign  in  the  coupling 
terms  since  the  restoring  forces  are  in  opposite  directions. 
If  the  distance  d  is  positive  which  means  the  level  of  water 
in  the  tank  is  higher  than  the  rolling  center,  the  magnitude 
gain  is  less  as  shown  in  Figure  5.3,  so  for  convenience  set 
d  equal  zero. 

Generally,  the  period  of  rolling  of  the  tank  (or 
natural  frequency)  must  be  designed  nearly  equal  to  the 
period  of  the  rolling  of  the  ship.   The  simplified  equations 

become : 

2  2  2 

$    +    as<}>   +    wn    <j>    -    y    a>n    6    =    wn    ty 

8  +  a.6  +  0)6  -  a).  <j>   =  0 


where 


a.  =  2t  u>_ 
s     s  s 


at  =  2ctwt 


0)   =  0)o  =  0). 

n    s    t 


By  manipulation  and  transform,  the  system  trans- 
mission function  of  the  rolling  angle  due  to  the  wave  slope 
force  is: 
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f(») 


wn(s2+atS+a)n) 


S4+(at+as)S3+(2u)^+atas)S2+W^(at+as)S+a)J(l-y) 


The   solution   of   the    equation    can   be   expressed    as 

an   amplification    factor  A(w)    and   e    the   phase    angle    lag   of 

the   forcing   function. 

2 


A   (<»)    = 


w. 


/x2   +    Y1 


1    Y 
X 


where      X 


••?   tan 

=     (wn    -    fa)    ) 


4/    2  2x 

li     Wp(Wn     -     03     ) 

7~2      272~      2    2 

(oj     -0)     )        +    OUfa) 


ua.w 
t   n 


-] 


Y        =   u[a:     +  —2      2~2+   2  T 

These  solutions  can  be  used  to  investigate  the 
effect  of  variation  of  the  tank  parameters  u.  and  y.  Calcu- 
lations can  be  carried  out  using  variation  ©f  each  parameter 
in  turn,  keeping  the  other  constant.  Adjust  both  parameters 
to  obtain  the  minimum  magnitude  gain.  These  values  of  para- 
meters   are    used    to    design   the    tank    stabilizer. 

The   computer   simulation   program   used   is    due    to 
Glavis    [7].      The   data    of   the    ship    in    the    previous    section 
were   substituted   into   the    equation   of   ship   and    tank   stabilizer, 
two   unknown    parameters    inthe   equation    are    ce»    and   y   as    variables 

Recal  1 


4(s)    =   T(s)    = 


o)2(S2+atS+.2) 


*%  S4+(a.+aJS3+(2aj2  +  a.ac)S2+a32(a.+aJS  +  u2(l-y) 


n      t   s 
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Numbering    the    data   of   the    ship   and   parameters    the 

equation   becomes: 

0.53023S2+0. 53023a. S+0. 28115 

T(s)  =  —. 3 S 2 

S4+( 0.046 5+at)S°+( 1.06047+0. 046 5at)S^+( 0.02466+ 

0.53023at)S+0.28115(l-y) 

For   the    sinusoidal    signal    the   equation   becomes: 

0.  281 15-0.  5302  3w2+j  0.5302  3a.  to 

T(jai)  =  -j 2 " 

a)  -(1  .  06047+0. 04651at)ui  +0.2811  5(1 -y)+j[(  0.024662 

+0.53023at)u-(0.046  5+at)u3] 

Then  9  «  ~ 

,  (0. 281 15-0. 53023w  r+ (0.5302 3a. w)  * 

[w4-(l  .06047+0. 046  51  at)a32+0.  281 15(l-y)]2+[(0.024662 
+0.53023at)u-(0.046  5+at)a)3] 

Since  the  size  of  the  tank  depends  on  the  free 
surface  coefficient  (y)  ,  if  y  is  large,  the  tank  will  be  big 
and  the  stability  of  the  ship  will  be  decreased  because  the 
metacentric  height  decreases.   So  the  magnitude  gain  of 
rolling  in  the  region  of  the  natural  frequency  decreases  but 
at  low  frequency  the  magnitude  gain  increases. 

Consider  the  computer  results.   Figure  5.4  shows 
curves  of  the  tank  damping  coefficient  versus  frequency,  with 
magnitude  as  curve  index  and  the  free  surface  coefficient  is 
constant.   On  Figure  5.4  the  free  surface  coefficient  is  0.05 
If  a  (tank  damping  coefficient  factor)  is  approximately  0.25 
the  height  of  the  resonance  peak  is  13  db ,  which  is  the 
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smallest  value  obtainable  for  these  conditions.   This  occurs 
at  a  frequency  of  about  0.72  rad/sec.   At  a  frequency  of 
about  0.18  rad/sec,  the  magnitude  is  only  1  db ,  and  is  even 
smaller  at  low  frequency. 

On  Figure  5.5,  the  free  surface  coefficient  is  0.1, 
the  line  of  constant  tank  damping  coefficient  factor  is  0.35 
to  give  a  resonance  peak  of  about  11  db  which  is  the  minimum 
obtainable.   Note,  however,  that  at  lower  frequencies  the 
magnitude  is  higher  than  for  the  condition  of  Figure  5.4. 

If  the  free  surface  is  set  equal  to  0.2  as  on 
Figure  5.6,  the  tank  damping  coefficient  should  be  0.45, 
which  gives  9  db  for  the  resonance  peak  and  2  db  at  low 
frequency  (0.5  rad/sec).   Thus,  if  the  free  surface  coeffi- 
cient is  increased,  the  magnitude  in  the  region  of  low 
frequency  also  is  increased,  though  the  height  of  the  reso- 
nance peak  is  decreased. 

The  curves  of  Figure  5.7,  Figure  5.8  and  Figure  5.9 
give  tank  damping  coefficient  factor  versus  magnitude. 
Frequency  is  the  curve  index,  and  the  free  surface  coefficients 
are  0.05,  0.1  and  0.2  respectively.   These  Figures  correspond 
to  Figure  5.4,  Figure  5.5  and  Figure  5.6.   By  inspection  they 
predict  the  same  optimum  values  for  the  damping  coefficient 
factor. 

Figure  5.10,  Figure  5.11  and  Figure  5.12  give  the 
Bode  diagrams  (plot  of  the  magnitude  versus  frequency)  with 
the  tank  damping  coefficient  factor  as  curve  index  and  the 
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free  surface  coefficient  factor  is  constant.   In  this  case 
the  Bode  diagram  cannot  show  the  best  value  of  the  tank 
damping  coefficient  factor  unless  that  value  has  been  plotted. 
However  the  Bode  diagram  is  good  for  consideration  of  the 
values  of  the  free  surface  coefficient  for  which  the  curves 
show  different  magnitudes  at  all  frequencies. 

By  inspection  of  the  curve  families  the  best  values 
of  the  tank  damping  coefficient  factor  and  the  free  surface 
coefficient  factor  are: 

the  tank  damping  coefficient  factor   =  0.25 
the  free  surface  coefficient  factor   =  0.05 
and/or 

the  tank  damping  coefficient  factor   =  0.35 
the  free  surface  coefficient  factor   =  0.1 
Note  that  the  magnitude  gain  value  for  low  fre- 
quencies increases  with  increases  of  the  free  surface  (y), 
which  means  increasing  tank  size.   This  results  from  the 
losses  in  roll  stiffness  due  to  the  tank  free  surface  effect. 
It  is  also  obvious  that  the  roll  response  at  low  frequencies 
increases  as  tank  size  increases. 

For  the  free  surface  tank  (as  reference  [8]),  the 
relation  between  the  parameters  and  the  tank  size  is: 

1  b3 

y    "   2AGM 

where  1   =   the  length  of  the  tank. 

b  =   the  width  of  the  tank  can  be  equal  or  less  than 

the  beam  of  the  ship. 
A  =  the  displacement  of  the  ship  in  tons. 
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GM 


metacentric   height 


The   other   suggestion   equation   is 


to 


h      = 


2     2 

t  b 

2  „ 

tt  g 


where     h   =  the  depth  of  the  water  in  tank. 

If  the  width  of  the  tank  is  assumed  to  be  equal  to 
the  beam  of  ship,  the  depth  of  the  water  in  the  tank  is  equal 
to  0.5  meter.   If  the  selected  free  surface  coefficient  factor 
is  0.1. 

1   =  0.82  meter. 

If  the  free  surface  coefficient  factor  is  equal  to 
0.05,  the  length  of  the  tank  is  reduced  to: 

1   =  0.41  meter 
which  can  save  space  in  the  ship.   And  the  other  parameters 
such  as  J.,  B.  ,  and  K.  can  be  calculated. 

The  time  response  to  an  initial  condition  is  shown 
in  Figure  5.13,  for  y  =  0.0,  0.05,  0.1,  and  0.2,  and  at  =  0.0, 
0.25,  0.35  and  0.45  respectively.   The  curve  for  y  =  0.0, 
without  antiroll  tank  stabilizer  is  very  oscillatory,  and  is 
damped  only  by  the  damping  of  the  ship  itself.   The  other 
curves  include  the  tank  stabilizer,  so  the  damping  is  higher, 
'and  the  curves  reach  steady  state  sooner.   The  time  constant 
of  the  undamped  system  is  about  42.5  seconds.   With  stabilizer 
added  the  setting  time  can  be  halved. 


58 


Figure  5.14  shows  the  time  response  to  a  sinusoidal 
input.   The  parameters  have  been  varied  as  for  Figure  5.13. 
The  frequency  of  the  input  is  set  at  about  the  natural  fre- 
quency (0.728  rad/sec).   For  the  ship  without  tank  stabilizer 
the  magnification  is  seen  to  be  excessively  large.   When  the 
tanks  are  added  the  system  response  is  seen  to  be  greatly 
reduced  in  magnitude. 

3.    U-Tube  Tank  Design 

The   U-Tude   tank   as    designed   in    this    chapter   is 
assumed    to   be   a    rectangular   tank.      The   mutual    inertia    term, 
J    .  ,   which    has    been    set    to   zero    in   preceding   sections,    is    not 
necessarily   zero.      The    value   of   J    .    depends    on    the    accelera- 
tion  of   the   water   relative    to    the    tank,    and    this    depends    on 
the    location   of   the    tank.      If   the    cross-duct   passes    below 
the   center   of   rotation,    J    .    is    a    positive    number  which    becomes 
smaller   as    the    tank    system    is    placed    higher,    as    shown    in 
Figure    5.15a    and   b. 

Consider   that   the   two    parameters    J.    and   J    .    are 


given   by 


Jt   =    pA0    R       J  A0    ds 
0   AT 


Jst  -  pAo  r2    5  h  ds 

u    R 
where    p      =   mass    density   of   the    fluid   in   tank. 

AQ   =    cross-sectional    area    of   vertical    leg    of   the    tank 
R     =    lever   arm   of   the    tank. 
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D  =  perpendicular  distance  from  "virtual"  center  of 
s   r 

rotation  to  the  tangent  at  any  point  on  the 

trajectary . 
A  =  the  cross-sectional  area  at  any  point  S. 
These  two  equations  can  be  related  as: 


J*t  -  fS  D*  A*  d5 


=   2h  + 


2  Dc  Ac 
s  s 


and  h  =  the  height  of  the  water  in  the  tank. 

Note  that  D  can  be  negative  if  the  tank  is  placed 

above  the  center  of  rotation.   So  J  +.  can  be  zero  if: 

st 


h  A, 


-D. 


Recall 

Jsi  +  Bs*  +  Ks<f>  +  Jste  +  Kste  =  kss^ 

•  •  •• 

ost*  +  Kst<f)  +  j  e  +  Bte  +  K.e  =  o 


By  using  the  definition  as  before  the  transformed 
equations  become. 

(S2+asS+o)2)(j,  +  (-  -jii  S2-  j^)e  =  ^ 

s       s 


(-  j5l  S2-u)2)<j>  +  (S2+atS+u)2) 
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K.         2 
Define  j^  =  3  (=yw  as  in  the  free  surface  tank) 

■^   =  X  =  2h  +  — ~ — S. 
dt  M0 


a4 


=  a 


then 


st  _ 


3X 
2 
n 


Using  Cramer's  rule,  the  ship  roll  transfer  function 
due  to  the  wave  slope  force  is: 


$(*> 


wn(S  +atS+o)n) 


( 1  -^-)  S4+  (  as+at )  S  3+  (  2wjj+atas  -  2X3 ) S  2+w2 (  at+as  )  S+w4  3a)2 


0) 


Inserting  the  values  of  the  parameters  of  the  ship, 
the  equation  becomes: 


l(s)  =  0.5302  3Sfc+0.5302  3aS+0.28n5 

*      [(1-1. 886X23)S4+(0.0465+a)S3+(l. 06047+0. 0465-2X3 )S2+ 

( 0.02466+0. 5 3023a) S+0. 2811 5-0. 530233 


then 
M2  = 


(0.28115-0.530  2  3a32)2+(0.5  30  23aco)2 

[(l-1.886X26)a34-(l  .06047+0.  0465 -2X3  )w2+0.  281 15-0.  530233]  2 

+  [(0.02466  2+0.5  3023aM0.0465+a)a)3]2 


Figure  5.16a,  b  and  c  are  plots  for  X=3  and  3=0.05, 
there  are  two  resonance  peaks  in  the  magnitude  curve,  one 
for  the  ship  itself  and  the  other  for  the  tank.   The  value  of 
a  at  the  minimum  peak  magnitude  is  about  0.45  and  only  one 
peak  magnitude  is  obtained.   The  point  P  in  Figure  5.16c 
corresponds  to  the  line  of  -2  db  in  Figure  5.16a. 
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If  the  tank  is  placed  as  high  as  possible,  the  value 

of  X  may  be  zero  or  negative.   If  X=0,  the  result  is  the  same 

2 
as  for  the  free  surface  tank  with  g=uw  '".   The  plot  of  X  =  -1.0 

and  3=0.05  is  shown  as  Figure  5.17a,  b,  and  c.   The  points  P 

and  Q  on  Figure  5.17c  are  the  intersection  of  the  a  lines, 

these  points  correspond  to  the  lines  of  the  frequencies  0.6 

and  0.9  rad/sec  on  Figure  5.17a.   In  this  frequency  region, 

if  the  value  of  a  is  decreased,  the  magnitude  also  decreases. 

Figure  5.17b  shows  that  minimum  peak  value  of  magnitude  is 

obtained  with  a=0.7. 

The  U-Tube  tank  design  equation  is: 

Kt  =  2PgAQ  R2 

where   P   =   mass    density   of   the   water. 

g   =   gravity   of   the   earth. 

Aq   =   cross-sectional    area    of   each    vertical    leg    of   the 

tank. 

R  =  distance  from  ship  center  line  to  center  of  a  tank 

vertical  leg. 

The  product  of  p  and  g  is  unity  for  the  water,  so  both 

can  be  eliminated  from  the  equation. 

And  the  cross-sectional  area  of  the  channel  between 

tanks  is  given  by: 

An  Rw2 
«   _   0   s 


s   n        h    2 
g  -  hws 
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Figure  5.3 


Frequency  Parameter  Plane  (Bode  Diagram).   Plot 
of  the  effect  of  the  water  level  in  Passive 
free-surface  tank. 
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Figure  5.4   Frequency  Parameter  Plane  for  Free^surface  Tank. 
Plot  of  Tank  Damping  Factor  vs  Frequency  with 
Free-surface,  y=0.05  and  M  as  curve  index. 
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Figure  5.5   Frequency  Parameter  Plane  for  Free-surface  Tank 
Plot  of  Tank  Damping  Factor  vs  Frequency  with 
Free-surface,  y=0.1  and  M  as  curve  index. 
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Figure  5.6   Frequency  Parameter  Plane  for  Free^surface  Tank 
Plot  of  Tank  Damping  Factor  vs  Frequency  with 
Free-surface,  y=0 . 2  and  M  as  curve  index. 
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Figure  5.7   Frequency  Parameter  Plane  for  Free-surface  Tank 
Plot  of  Tank  Damping  Factor  vs  Magnitude  with 
Free-surface,  y=0.05  and  to  as  curve  index. 
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Figure  5.8 


Frequency  Parameter  Plane  for  Free-surface  Tank 
Plot  of  Tank  Damping  Factor  vs  Magnitude  with 
Free-surface,  y  =  0.1  and  w  as  curve  index. 
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Figure  5.9 


Frequency  Parameter  Plane  for  Free-surface  Tank 
Plot  of  Tank  Damping  Factor  vs  Magnitude  with 
Free  Surface,  u  =  0.2  and  w  as  curve  index. 
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Figure  5.10   Frequency  Parameter  Plane  for  Free-surface  Tank, 
Bode  Diagram,  Plot  of  Tank  Damping  Factor  as 
curve  Index  with  Free-Surface,  y=0.05. 
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Figure  5.11 


Frequency  Parameter  Plane  for  Free-surface  Tank 
Bode  Diagram,  Plot  of  Tank  Damping  Factor  as 
Curve  Index  and  Free-surface,  u=0.1. 
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Figure  5.12 


Frequency  Parameter  Plane  for  Free-surface  Tank 
Bode  Diagram,  Plot  of  Tank  Damping  Factor  as 
Curve  Index  and  Free-surface,  y=0.2. 
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Figure    5.13      Initial    Condition    Response   with    Different    Value 
of   Tank    Damping    Factor   and    Free-surface. 
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(a) 


(b) 
Figure    5.15      U-Tube   Tank 
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Figure  5.16a 


Frequency  Parameter  Plane  for  U-Tube  Tank. 
Plot  of  Tank  Damping  Factor  vs  Frequency 
with  X  =  3,  3=0.05  and  M  as  Curve  Index. 


76 


Figure  5.16b   Frequency  Parameter  Plane  for  U-Tube  Tank 
Plot  of  Tank  Damping  Factor  vs  Magnitude 
with  X  =  3,  B  =  0.05  and  co  as  Curve  Index. 
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Figure  5.16c 


Frequency  Parameter  Plane  for  U-Tube  Tank 
Bode  Diagram,  Plot  of  Tank  Damping  Factor 
as  Curve  Index  with  X  =  3  and  3=0.05. 
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Figure  5.17a 


Frequency  Parameter  Plane  for  U-Tube  Tank 
Plot  of  Tank  Damping  Factor  vs  Frequency 


with  X 


1,  3=0.05  and  M  as  curve  index 
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Figure  5.17b   Frequency  Parameter  Plane  for  U-Tube  Tank 
Plot  of  Tank  Damping  Factor  vs  Magnitude 
with  X  =  -1 ,  6=0.05  and  u  as  Curve  Index. 
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Figure  5.17c 


Frequency  Parameter  Plane  for  U-Tube  Tank, 
Bode  Diagram,  Plot  of  Tank  Damping  Factor 
as  Curve  Index  with  X  =  -1  and  3=0.05. 
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This  is  for  the  rectangular  U-Tube  tanks.   With  the 
equations  as  above,  the  values  of  cross-sectional  area  of 
the  vertical  leg  of  the  tank,  AQ,  the  cross-sectional  area 
of  the  channel,  A  ,  can  be  calculated  and  tank  designed. 

The  depth  of  the  water  in  tank,  h,  and  the  distance 
between  the  rolling  center  and  the  channel  of  the  tank,  D  , 
can  be  calculated  by: 


X  =  2  h  + 


2  Ds  As 

An 


In  the  case  of  X  =  3,  if  the  tank  is  located  below 
rolling  center,  h  should  be  less  than  1.5. 

E.g.,  if  X  =  3,  3=0.05,  the  tank  damping  factor  is 
0.45.   Then  the  tank  size  is: 

2    Kt   ^Js         4 
A0  R  =  2P?  =  2P?  =  32-25  m 

if  R  =  4,  the  width  of  the  tank  is  about  1  meter.   Therefore 

AQ  =  2.013  m2 

2 
and  if  h  =  1.3  meters,  the  result  of  A  is  0.4686  meter 

and  Ds  =  0.859  meter. 

E.    ACTIVE  TANK  STABILIZATION 
1 .   Analysis 

Acti vated-anti rol 1  tanks  are  considered  in  detail. 
Particular  attention  is  paid  to  the  design  of  the  tank  control 
system.   The  equations  of  motion  are  the  same  as  those  of  the 
ship  with  passive  tank  stabilizer,  but  must  be  modified  to 
include  the  pump  used  in  the  activated  tank.   Active  tank 
control  can  best  be  selected  on  the  basis  of  minimizing  the 
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response  of  the  ship  to  an  impulsive  roll  velocity,  and  it 
will  be  shown  that  the  active  tank  is  significantly  better 
in  reducing  the  roll  of  ship  than  the  ordinary  passive  antiroll 
tank.   The  active  tank  stabilizer  studied  uses  the  U-Tube  tank 
design. 

Active  tank  analysis  parallels  that  of  passive  tanks 
and/or  the  double  pendulum.   The  tank  system  is  activated  by 
inserting  a  force  (such  as  a  pump  system)  into  the  second 
equation  of  passive  tank  system.   The  equation  can  then  be 
written  as: 

ship:   Js<|»  +  BS<J>  +  Ks4>  +  Jst'e  +  Kst0  =  K^ 

tank:   J  t<|>  +  K  t<|>  +  Jte  +  Bte  +  Kt8  =  F 

where  Fp  =  the  force  due  to  pump  system. 

The  effect  of  the  roll  sensing  control  on  the  force 

produced  by  the  pump  system  can  be  expressed  as: 

2 


Fp(s)  - 


KpfG^GgS  +  GgS'H 


1+  xS 


where  K_  =  the  gain  of  the  control  pump. 

G,  =  roll  angle  sensitivity. 

Gp  =  roll  rate  sensitivity. 

G_  =  roll  acceleration  sensitivity. 

t   =  time  delay  constant. 
For  convenience,  we  set  x=0  in  the  computation,  since 
the  time  delay  is  small. 
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The  active  tank  system  can  be  shown  in  block  diagram 


form  as 


PASSIVE 


* 


F^s) 


^9 


* 


PUMP 


F2(s) 


CONTROL 


where  F^(s)  = 


F2(s) 


F0(s)  = 


transfer  function  of  the  ship  rolling  due  to 
wave  slope  alone  which  is  the  same  as  for  the 
passive  tank  system. 

transfer  function  of  the  ship  rolling  due  to 
pump  alone,  which  can  be  expressed  as: 

3   


4  „  2 


(l-y-)S4+(as  +  at)S3+(2o)^+atas-2X3)S2+a3;(at+as)S+o)n-Ba)n 


2.   System  Design 

From  the  analysis  in  this  section,  the  system  per- 
formance can  be  studied  using  computer  simulation.   The  pump 

■    .    K 

coefficient  factor  K   (K  =  -*£■)  is  one  variable  parameter 

p   p   Jt 

which  may  be  called  a.   It  is  varied  as  the  pitch  ratio  of 
the  propeller  of  the  pump,  and  can  be  assumed  to  be  varied 
from  zero  to  one  and  zero  to  minus  one.   For  convenience  the 
actual  system  may  be  designed  with  only  positive  variation  and 
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the  direction  of  the  water  motion  between  tanks  is  controlled 
by  valves. 

From  the  block  diagram  the  system  transmission 
function  is: 


f<*>  = 


o)2[S2+at    S+co2] 


(l-y-)S4+(at+a    )S3+(2aJ^+ata    -2XB-K'G3B)S2+(W2at+^ae 


where 

V 

Is) 

=    Kp(G 

Define 

X 

=    3 

3 

=    0.05 

a+ 

=    0.4 

and 


K. 


=    a 


These  values  for  a  passive  tank  gave  a  peak  resonance 
of  minimum  magnitude  at  about  10  db  at  0.6  rad/sec  and  is  1  db 
at  the  lower  frequencies.   Using  these  values  the  equation 
becomes : 


(s)  = 


0.530  2  3S^+0.21209  2S+0.28115 


0. 1 51 3S4+0.4465S3+(0. 7791 -0.05aG3)S2+(0. 237-0. 05aG2)S 


+0. 25464-0. 05aG 


1 


Figure  5.18  shows  the  Bode  diagram  for  G,=l, 
G2=G3=0.   If  a  is  more  positive,  the  frequency  and  the  magni 
tude  of  the  resonance  peak  are  slightly  decreased.   At  the 
low  frequencies,  the  magnitude  is  increased.   If  a  is  more 
negative,  the  result  is  reversed.   To  obtain  the  minimum 


85 


Figure  5.18 


Frequency  Parameter  Plane  for  Activated  Tank. 
Bode  Diagram,  Plot  of  Pump  Coefficient  Factor 
as  Curve  Index  (G,=l,  G9=GQ=0)  and  X=3. 


'1 


2  u3 
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magnitude  gain,  at  the  lower  frequencies  from  0  to  0.6  rad/sec 
a  is  -1,  for  the  higher  frequencies  a  must  be  +1 .   So,  the 
pitch  ratio  of  the  propeller  of  the  pump  may  be  kep  constant 
and  the  direction  of  the  water  motion  reversed  by  switching 
the  valves.   Then  the  maximum  magnitude  is  7  db  at  0.6  rad/sec 
and  -1  db  at  the  lower  frequencies. 

The  characteristics  of  the  active  device  (pump, 
valves,  sensor)  can  be  treated  as  the  variation  of  the  value 
of  the  parameter  a  versus  frequency.   This  value  is  required 
to  switch  at  the  resonant  frequency  as  shown  in  Figure  5.19. 


1 

a 

0 

0.6 

W  rad 

-1 

-1 

Figure  5.19  The  variation  of  control  device. 
(X  =  3,  G^l  ,  G2  =  G3  =  0) 

For  G-.  =0 . 5 ,  Gp  =  l  .25  and  G3  =  2.4  the  results  are  shown 
as  Figure  5.20a,  b  and  c.   Figure  5.20a  is  a  plot  of  a  versus 
frequency,  and  magnitude  is  the  curve  index.   When  a  is  more 
negative  the  magnitude  is  seen  to  be  smaller.   If  a=l  the 
maximum  peak  magnitude  is  obtained  at  the  frequency  of  about 
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Figure    5.20a 


Frequency  Parameter  Plane  for  Activated  Tank 
Plot  of  Pump  Coefficient  Factor  vs  Frequency 
with  M  as  Curve  Index  (G,=0.5,  G«=1.25, 
G~=2.4)  and  X=3.  c 
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Figure  5.20b   Frequency  Parameter  Plane  for  Activated  Tank 
Plot  of  Pump  Coefficient  Factor  vs  Magnitude 
with  u  as  Curve  Index  (G,=0.5,  G0=1.25, 
G3=2.4)  and  X=3.        '       L 
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Figure  5.20c 


Frequency  Parameter  Plane  for  Activated  Tank. 
Bode  Diagram,  Plot  of  Pump  Coefficient  Factor 
as  Curve  Index  (G,=0.5,  G9=1.25,  G~=2.4)  and 
X=3.  '       c  6 
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0.6  rad/sec.   These  results  correspond  to  Figure  5.20b  which 
shows  the  frequency  curve.   At  the  same  frequency  of  0.6 
rad/sec,  both  figures  show  maximum  magnitude  for  all  values 
of  a  from  -1  to  +1.   Figure  5.20c  is  the  Bode  diagram  and 
shows  that  there  is  no  phase  lag  between  ship  and  active  tank, 
and  if  a  is  more  negative  the  magnitude  is  smaller.   So,  at 
a=-l  ,  it  produces  the  smallest  magnitude  at  all  frequencies, 
with  the  peak  equal  5  db  at  0.6  rad/sec. 

This  result  seems  to  be  the  best  since  all  the 
parameters  of  the  pump  system  (pitch,  sensing  gain  control, 
and  valves)  are  fixed. 

From  both  results  as  mentioned  above,  the  system 
can  be  designed  as  a  linear  device. 

If  there  is  no  limit  on  the  height  at  which  the 
tank  can  be  located,  it  is  possible  to  choose  the  height,  so 
that  X=0.   If  this  is  done,  by  inspection  the  value  of  3=0.05 
and  a. =0.3  are  obtained.   Substitute  these  values  into  the 
equation  of  the  system  which  becomes: 


(s)  = 


0.5302 35^+0.159  0695+0. 2811 5 


S4+0.3465S3+(1 .07436-0. 05K  G- ( S2+( 0 . 18373-0 . 05K  G,)S 

p    3  p    l 


+0. 25464-0. 05K    G] 

Figure    5.21a    and   b    show   the    results   with    G-,=l    and 
G2=G3=0.      Figure    5.21a    is    a    plot   with   magnitude   as    curve 
index;    point   Q    is    the    maximum   magnitude   with    a=-0.65    at    a 
frequency    of   0.728    rad/sec.      If   a=0.15,    the    maximum  magnitude 
of   9    db    is    obtained   which    is    the    minimum    peak.      Figure    5.21b 
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Figure  5.21a   Frequency  Parameter  Plane  for  Activated  Tank. 
Plot  of  Pump  Coefficient  Factor  vs  Frequency 
with  M  as  Curve  Index  (G-,=l,  G2  =  G'3=0)  and  X  =  0 
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Figure  5.21b 


Frequency  Parameter  Plane  for  Activated  Tank. 
Bode  Diagram,  Plot  of  Pump  Coefficient  Factor 


as  Curve  Index  (G-j  =1  ,  G 


:G3=0) 


and  X=0 
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shows  the  Bode  diagram.   To  obtain  the  minimum  magnitude,  in 
the  frequency  region  between  0.6  and  0.85  rad/sec  the  value 
of  a  must  be  maintained  at  +  1  and  for  the  other  frequencies, 
at  -1.   There  are  two  magnitude  peaks  of  approximately  7  db 
at  the  intersection  of  a=+l  and  a=-l .   For  the  lower  frequen- 
cies the  magnitude  is  -1  db. 

From  the  above  result  the  active  device  parameter 
value,  a,  can  be  plotted  versus  frequency  as  in  Figure  5.22. 
This  figure  shows  two  values  of  a  at  the  frequencies  between 
0.6  and  0.85,  the  straight  line  corresponds  to  the  minimum 
magnitude,  and  the  curved  line  corresponds  to  a  constant 
magnitude  in  that  region. 


a)  rad/sec 


Figure  5.22  The  Variation  of  Active  Device 
(X  =  0,  G^l  ,  G2  =  G3  =  0) 

Figure  5.23  shows  the  Bode  diagram  of  the  system 
for  G,=0.5,  G2= 1 . 25 ,  and  G3=2.4.   For  minimum  magnitude  with 
the  control  device,  the  required  variation  in  a  is  shown  in 
Figure  5.24. 
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Figure  5.23  Frequency  Parameter  Plane  for  Activated  Tank. 
Bode  Diagram,  Plot  of  Pump  Coefficient  Factor 


as  Curve 
X  =  0. 


Index  (G.j=0.5,  G2  =  l  .25, 


G3=2 


4)  and 


95 


a)  rad/sec 


Figure  5.24  The  Variation  of  Control  Device 
(X=0,  G^O.5,  G2  =  l  .  25  G3  =  2.4) 

Figure  5.25,  which  is  the  Bode  diagram,  compares 
the  results  for  an  undamped  ship  and  for  a  ship  damped  by  a 
tank.   All  have  only  one  resonance  peak.   The  magnitude  of 
roll  in  the  undamped  ship  at  resonance  is  the  highest,  damping 
clearly  reduces  the  amplitude.   However,  in  the  passive  tank 
system  the  magnitude  at  lower  frequencies  is  slightly  higher, 
and  the  curves  show  that  the  U-Tube  tank  is  more  advantageous 
than  the  free-surface  tank.   The  active  tank  system  gave  the 
most  advantageous  behavior  as  shown. 
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Figure    5.25      Compare   the    Results    in   Magnitude    Gain 
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VI.   CONCLUSION 

Analysis  techniques  can  be  employed  effectively  in  system 
design  particularly  when  an  exact  analytical  solution  is  not 
available  and  where  a  trial  and  error  method  is  the  only 
solution  possible. 

The  frequency  parameter  plane  techniques  have  been  used 
in  this  thesis  to  analyze  and  design  a  system.   The  system 
elements  which  may  be  linear  or  nonlinear  were  represented 
by  adjustable  parameters.   To  select  a  stable  system,  the  root 
parameter  plane  was  used  and  further  mapped  onto  the  frequency 
parameter  plane.   The  goal  was  to  obtain  the  best  set  of  para- 
meters for  a  stable  system. 

Rather  than  adjust  parameters  for  every  trial,  one  graph 
was  made  and  the  parameter  value  corresponding  to  the  best 
curve  was  chosen.   This  appreciably  reduced  the  number  of 
trials  and  the  computation  time  in  the  overall  design. 

In  the  system  which  is  presented  the  tank  stabilizer 
design  might  not  be  optimal  ,  however  further  effort  on  para- 
meter adjustments  by  iteration  of  the  described  technique 
would  improve  the  design. 

The  tank  stabilizers  of  a  ship  were  designed  in  this 
thesis  by  using  the  frequency  parameter  plane  technique,  they 
were  wery   successful  and  effective.   If  the  relationship 
between  the  designed  parameters  and  tank  size  are  true,  the 
result  of  designing  by  using  the  method  would  be  feasible  and 
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believable.   In  passive  tank  design,  the  tank  size  depends 
on  the  free  surface  factor.   The  tank  damping  factor  depends 
on  the  skin  of  the  material  of  the  tank,  which  would  be  tested 
to  obtain  its  value  from  the  analysis.   The  value  of  X  in 
U-Tube  tank  would  have  to  be  chosen  within  the  range  which 
maintains  ship  stability. 

In  the  activated  tank  design,  the  gain  control  constants 
were  varied  by  trial  and  error  to  adjust  the  damping.   Phys- 
ically only  one  parameter  can  be  varied,  the  pitch  radio  of 
the  propeller  of  the  pump. 
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